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Abstract: 

We compute the effects of a stochastic background of gravitational waves on multiply imaged 
systems or on weak lensing. There are two possible observable effects, a static relative deflec- 
tion of images or shear, and an induced time dependent shift or proper motion. We evaluate 
the rms magnitude of these effects for a COBE normalized, scale-invariant spectrum, which 
is an upper limit on spectra produced by inflation. Previous work has shown that large-scale 
structure may cause a relative deflection large enough to affect observations, but we find that 
the corresponding effect of gravity waves is smaller by ~ 10 4 and so cannot be observed. 
This results from the oscillation in time as well as the redshifting of the amplitude of gravity 
waves. We estimate the magnitude of the proper motion induced by deflection of light due to 
large-scale structure, and find it to be ~ 10 -8 arcsec per year. This corresponds to ~ 50 km/s 
at cosmological distances, which is quite small compared to typical peculiar velocities. The 
COBE normalized gravity wave spectrum produces motions smaller still by ~ 10 2 . We con- 
clude that light deflection due to these cosmological perturbations cannot produce observable 
proper motions of lensed images. On the other hand, there are only a few known observational 
limits on a stochastic background of gravity waves at shorter, astrophysical wavelengths. We 
calculate the expected magnitudes of the effects of lensing by gravity waves of such wave- 
lengths, and find that they are too small to yield interesting limits on the energy density of 
gravity waves. 



I Introduction 

Events in the early Universe may have left a stochastic background of gravitational waves (GW's). In 
particular, a generic prediction of inflation is a relic spectrum of GW's Detecting these elusive 
remnants would not only establish this prediction of general relativity, but also serve as a critical test for 
inflation. While the predicted background may be too weak for direct detection ||, it could be detected 
indirectly through its effect on light propagation in the Universe. Even if the effects of GW's cannot be 
distinguished observationally from other effects, observers who assume no GW's might reach incorrect 
conclusions about the distribution of matter in the Universe. 

Gravitational lensing is one of the most promising methods of mapping the distribution of matter at 
cosmological distances. Detailed observations of multiple images of quasars have been used to try to 

Electronic address: barkana@arcturus.mit.edu 



1 



reconstruct the lensing mass distribution (e.g., ||). It has also long been recognized that measurements 
of the time delay between images can be used to determine the Hubble constant [|J. Gravitational lenses 
and sources, however, typically lie at significant redshifts. Light rays are thus deflected by large-scale 
structure (LSS) and GW's as they traverse the cosmological distance to the observer, and these deflections 
may change the simple lensing picture. 

GW's may be produced by many sources. Astrophysical sources, such as close binary systems which 
include a neutron star or black hole, radiate GW's, and numerous individual sources may superpose to 
create a stochastic background. At the Planck time, quantum fluctuations in the metric are significant 
and may produce gravitons. Phase transitions in the universe may lead to topological defects such as 
cosmic strings, which generate GW's. A period of inflation may leave behind a significant amount of 
GW's. Whatever the source, any spectrum which extends over wavelengths comparable to the present 
horizon would contribute to the quadrupole anisotropy of the cosmic microwave background (CMB) |J. 
Such a spectrum is therefore limited by the anisotropy measured by the COBE DMR experiment ||. 
For our calculations we adopt a scale- invariant primordial spectrum, i.e., one which has constant energy 
density per logarithmic frequency, which we assume produces the entire measured quadrupole anisotropy. 
Inflationary models predict slightly tilted spectra which are responsible only for some fraction of the 
anisotropy ||, |2|, and so are generally weaker than our adopted case. 

In inflation, GW's are produced in conjunction with density fluctuations. The initial nearly scale- 
invariant power spectrum of density fluctuations evolves as modes reenter the horizon after inflation, 
and as structure later forms in a universe dominated by dark matter. The present spectrum is strongly 
constrained by galaxy and cluster surveys, and can be used to study the effects of LSS on lensing. The 
induced effects are small but potentially observable. In weak lensing, the effect is a coherent distortion of 
background galaxies by an ellipticity of the order of a few per cent M, M. In strong lensing, the primary 
effect is an external shear which may be significant for observed four-image systems g[fi}. 

In general, the influence at a given time of a weak metric perturbation on light propagation is simply 
described by two effects. Their magnitudes were estimated for LSS in Ref. |1C], which we summarize 
here. The first effect is a constant deflection, the same for all nearby light rays. This deflection simply 
displaces the "true" angular position of an observed lens or source, and is not directly observable. In 
the case of LSS, deflections from coherent structures of size ~ 1 Mpc add up in a random walk, giving 
an overall deflection of order a few arcminutes at redshift 1, which scales as the square root of comoving 
distance r. The second effect is a relative deflection between nearby light rays, which produces a focusing 
and shear with observable effects on weak and strong lensing. For two rays at initial angle 9, each 
coherent structure at a distance r causes a relative deflection proportional to their separation of ~ r6. 
The additional random walk gives a relative angular fluctuation of pa 0.07 6 at redshift 1, which scales as 

r 3/2_ 

It was suggested in Ref. JTI| that gravity waves could significantly affect the time delays in a multiply 



imaged system. It was later pointed out [I2| that a correct analysis must include the lensing constraint, 
i.e., the fact that image rays in the presence of GW's follow different paths than for no GW's, so that 
all rays go from the source to a common destination, the observer. These later authors also showed that 
both LSS and GW's have no observable effects on lensing, to lowest order. However, in their lowest order 
expansion they assumed that two image rays that are observed at an angular separation are separated 
by a distance of exactly r6 on the lens plane at a distance r. In other words, they neglected the relative 
deflection between light rays, and therefore only included an overall, constant deflection due to LSS or 
GW's. 
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We can easily see why this assumption leads to no observable effects. In the absence of metric pertur- 
bations, we can write the lens equation for a thin lens as (e.g., |L3|]) (3 = 6 — a\ ens (9), where 9 and (3 are 
the image and source angles, respectively, and c?i en s is the scaled deflection angle, which is determined 
by the mass distribution of the lens. If we neglect relative deflections, then LSS or GW's can only cause 
an angular shift a^ between the observer and the lens, and a shift a^s between the lens and the source. 
Then the lens equation becomes (3 = 9 — a\ ens (9) + as, where 9 is now measured relative to the observed 
(and shifted) lens position, and as involves «l and &ls (see Sec. Ill for the full details). The constant 
(i.e., ^-independent) deflection as has no effect on any observables of the lens system (e.g., Jll|), since (3 
is not directly observable . Fermat's principle then implies that the lens equation must be equivalent to 
dAt/d9 = at fixed (3, where At is the relative time delay. There is thus no observable effect on the time 
delay, either, since it can be derived from the lens equation, up to (unobservable) ^-independent terms. 

This approximation of neglecting the relative deflection may not be a good one. Indeed, such deflection 
can have observational consequences, which may be sufficiently large to detect in the case of LSS || [To| ] . 
In this paper, we compute the rms total and relative deflections between light rays induced by a scale- 
invariant stochastic background of GW's. Unlike LSS, GW's oscillate with time, and so the effect of 
short wavelength modes does not amplify, as light rays deflect one way in crests and the opposite way in 
troughs. In addition, the energy density and thus also the amplitude of sub-horizon GW's redshift away 
as the universe expands. The lensing effect is thus dominated by wavelengths on the scale of the distance 
to the source. Each such mode acts as a single coherent structure, and so both the total and relative 
deflections due to GW's scale approximately linearly with distance. The effect of different modes must 
be convolved with a particular power spectrum and include the above-mentioned decay of each mode as 
the universe expands. We find simple integral expressions for the scale-invariant spectrum. The total 
and relative deflections are smaller than those caused by large-scale structure by factors of the order of 
10 2 and 10 4 , respectively. We do not need to explicitly set up the lens equation, since the rms shear in 
the lens equation is directly related to the rms relative deflection of light rays, which we calculate. This 
fact was demonstrated for LSS in Ref. Q, and we give a general proof in Sec. Ill below. Our results 
imply that the static effects of the GW spectrum on lensing are negligible compared to those of LSS, and 
cannot be detected in practice. 

In addition to the static effects of LSS and GW's on lensing, it is possible that the fluctuation in the 
induced deflection with time would be directly manifested as an observed proper motion of images. In 
other words, the sources do not really move but the light rays from the sources are deflected and so 
the sources appear to move. We find that even LSS can only produce motions of order 10 -8 arcsec per 
year from this effect. This corresponds to ~ 50 km/s at a distance of a Gpc, and the effect of GW's is 
smaller still by a factor of ~ 10 2 . Since typical peculiar velocities are much larger, the proper motion 
induced by deflection of light due to LSS is unobservable, and the same is true for the COBE-normalized 
scale-invariant spectrum of GW's. 

However, we may try to use shear or proper motions of imaged sources to improve existing limits 
on stochastic GW's at a range of astrophysical wavelengths. There are only a few such limits known: 
Single-pulsar timing yields Q\ < 1 x 10 -8 at A ~ 2 p c |l4| , |l5| , binary pulsar timing implies Q\ < 0.04 
over A ~ 2 pc to 1 kpc and Q\ < 0.5 up to 10 kpc [14, [l6|], and the observed angular correlation function 
of galaxies sets a limit of Q,\ < 1CP 3 over A ~ 100 kpc-100 Mpc [17]. These limits apply to any stochastic 
background of GW's, whether cosmological in origin or generated at low redshift as a superposition of 
many discrete sources. For a cosmological spectrum that existed at early times, there are also big bang 
nucleosynthesis constraints of Q\ < 10~ 4 for A < 100 pc [18] and CMB limits of Q\ < 10~ 12 at Horizon 
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wavelengths (from COBE) and £1\ < 10 for A > 1 Mpc from small-scale anisotropy [jn 

In Ref. |l7j it was suggested that highly magnified lensed sources could increase the sensitivity to de- 
tecting proper motions due to GW's. The angular deviations induced by GW's produced by an individual 
source were discussed in Ref. [2C]. Ref. [^] considered detecting proper motions (of unlensed sources) 
due to GW's through VLBI measurements, but our approach is simpler than theirs. For an image of a 
lensed source, only an angular deflection of the source relative to the lens is easily observed, and we find 
that this relative motion is small when we assume an isotropic GW background. Thus we do not find an 
interesting limit on the energy density. 



II Formalism 



In this section we review the formalism describing gravity waves, their cosmological evolution, and their 
effect on lensing, as well as the usual formalism of gravitational lensing. We work in the framework of a 
flat Robertson- Walker metric with small-amplitude tensor metric fluctuations. For weak perturbations, 
we can consider the effect of GW's without including LSS, since the cross terms between them would be 
of higher order. In comoving coordinates we can write the line element as 



a 2 (r)[-dr 2 + (<% + h^dx'dx^ . 



(1) 



Here r is the conformal time, a(r) the expansion factor, and we have set c = 1. We expand the metric 
perturbation in plane waves {k = 2ir/\), 



h lm (x,r)= d 3 k h n (k,r)ef m (k)e 



-ik-x 



(2) 



where ef m is the polarization tensor which depends on the direction k (I and m are spatial indices ranging 
from 1 to 3, while n ranges over the polarization components +, x). For a wave propagating in the 
z-direction, the nonvanishing components are in the x-and y- rows and columns: 
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For other propagation directions k, we rotate e — > ReR T , with R the standard 3x3 rotation matrix. 

GW's with a given wavevector k are produced during inflation and then stretched outside the horizon. 
The amplitude is constant outside the horizon, but once a mode reenters its energy redshifts as a~ 4 . 
For the inflationary spectrum the effect of very short wavelength modes is negligible, and so we can 
assume that all modes enter during the matter-dominated era, for which the exact time evolution is 
given jl], H by a spherical Bessel function, 3ji(/cr)/(A;r). This time evolution is also correct for all modes 
long after matter-radiation equality. Inflation produces Gaussian, stochastic perturbations. The Fourier 
components have zero ensemble mean and a covariance 



< h i (k,T 1 )h j (q,T 2 ) >= A T k~ 



kr\ 



3ji(fcr 2 ) 
kro 



8 3 (k + q)8., 
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for the scale-invariant k 3 spectrum. Note that we do not assume the short-wavelength approximation 
h l (k,T) oc a _1 (r) e lkr . The contribution to O at the present (averaged over several periods) is 

° A = IHnT = T^ T(fcT0) ' (4) 

where ro = ZHq 1 is the present value of r, and throughout we set Hq = 75 km sec" 1 Mpc -1 . Normalization 
to the full CMB quadrupole anisotropy gives A-j = 6 x 10 -11 . 

Consider a photon emitted from a source toward an observer at the origin, with the photon's final 
direction defined as (minus) the z-axis. We use r to denote values of the z-coordinate (with zs denoting 
the source redshift, not its z-coordinate). GW's affect the distance-redshift relation, but this effect is 
separate from that of the angular deflections which we are interested in, and it introduces only small 



additional corrections in these quantities [17]. We can thus neglect this effect, and assume that the 
photon path obeys r(r) = tq — t. In a flat, matter-dominated universe, rs = 2Hq 1 [1 — (1 -Ms) -1 / 2 ]. The 
components perpendicular to the z-axis of the photon direction obey [fiT 



(r) ~ -r-(ro) = h zi (r) - h zl (r ) - - / V»MO dr' . (5) 



Integrating this we find, for the perpendicular components of the position (with respect to x(tq) = 0), 



At) 



~(t' - T)Vih zz (r') + h zi (T ) - h zi (r') 



dr' . (6) 



We define a (two-component) angle (3 l = x 1 (t) /t(t). 

In gravitational lensing with a primary thin lens at a distance tl (but no LSS or GW's) the lens 
equation is (e.g., fl3|l ) 

P = 0-a lena (6) , (7) 

where 8 is the observed image angle, (3 is the source angle (defined as xs/rs, in terms of the perpendicular 
position of the source), and c?i e ns is the deflection angle scaled by r^s/rs (we define r^s = rs — 7*l). In 
this case, the fiducial z-axis is defined to be in the observed direction of the lens. The distortion of the 
image of a small source is given by the inverse of the Jacobian matrix 

w = 6 "- ' (8) 

where is also termed the shear tensor of the lens. 



Ill Shear induced by GW's on lensing 



In this section we follow the approach used for LSS in Ref. [10|, i.e., we compute some of the same 
quantities for GW's and compare the results. As stated in Sec. I, we do not need to include a lens 
explicitly, as we now justify. In the presence of a metric perturbation, but without a primary lens, the 
lens equation has the form (3 = 6—aos{6)i where aos results from the accumulated deflection between the 
observer and the source. As defined in Sec. II, the shear tensor for an image at 9 due to the perturbation 
is i™ = da l Q S (6)/d6 : ' . On the other hand, the relative deflection at 9 between two rays separated by 
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a tiny angle 7 is aos(G + 7) — ^os{&)- We denote the rms of this quantity by (TAp- We average over 
directions of 7 (which for this calculation is equivalent to fixing 7 and assuming that F %3 is isotropic) and 
take 7 — > 0, obtaining the relation ^] 

(fM)'^^), (9) 

all evaluated at position 8. Thus, CTA/3/7 yields an estimate of the magnitude of the shear tensor. Indeed, 
it fully characterises rms values of F* 3 , since for an isotropic field 

(FijFu) = I [S i:j S kl + SikSji + S u S jk ] (F mn F mn ) . (10) 

o 

If we also include a primary lens, in the lens equation we simply add up all deflections linearly, 
assuming all deflections are small. For the primary lens alone, we have Eq. (|7|). Figure |l] shows this setup 
schematically. In the presence of a metric perturbation, we trace a light ray that is observed at r = to 
come from the direction 8, back to the source. We find a different form for the lens equation: 

= 0- a { of ~ 3 lew (9- agi) • (11) 

—(2 3) 

Here a ' s refers to the integrated deflection caused by LSS or GW's along the paths labeled 2 and 3 

n (2 3) (2} 

in Figure [I], defined so that the total induced change in x(ts) equals rsa 's ■ Similarly, rLa OL is the 
induced change in x(t£). In integrating the deflections along the unperturbed paths 2 and 3, we are 
assuming that the relative deflections due to LSS or GW's are small compared to 8 and aien S) which is 
true for the cases which we consider below. 

When the perturbation is included, 8 is no longer an observable, since it is measured with respect to 
the unperturbed position of the lens. The observed position of the lens (whose actual position has not 
changed) is now 8\ ens = Oq L , so the lens equation in terms of the observable 8' = 8 — 8\ ens is 

= 9' + o%\ - a%f - a lens (8' + agi - o$ L ) . (12) 

If we now calculate the shear tensor resulting from equation ( |l2[ ) , it will contain the shear of the primary 
lens, shear terms from the perturbation, and also cross terms. For simplicity, in the case of GW's we 
only estimate one characteristic magnitude, that of the shear resulting from a g , by evaluating the 
corresponding cja/3/7- In Ref. || all the different shear terms were studied for LSS, showing that 0A/3/7 
indeed estimates the relative magnitude of the various corrections due to LSS. Since we find that aA/3/7 is 
much smaller for GW's, we do not have any motive to explore equation (|l^) further in this case. Instead 
of the path (2,3), we may use a straight path from r = to r = r$ to evaluate the rms of various 
quantities in this section, since aq e ns <C 1 an d so the components of vectors and tensors as well as the 
relative distances of points along the path (both of which enter into the rms calculations) are unchanged 
(except for 0(a\ ens ) corrections). Thus we only need to consider the effect of GW's in the absence of a 
primary lens. 

Consider first a single light ray with 8 = 0. In the absence of GW's (or LSS) it would follow the 
straight line x 1 {t) = for all r. We now include the effect of GW's, and compute the rms fluctuation in 

2 Repeated indices are summed over the x and y directions. There is no distinction between upper and lower indices. 
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Figure 1: Sketch showing positions of the observer, lens and source, as well as an image ray and several 
comoving distances. 



i- - V 1 / 2 

the photon's perpendicular displacement at the source, up = ((3(ts) ■ /3(rg) ) ■ This is a measure of the 
common deflection of all image rays, and is therefore not observable, but it is useful for the calculations 
that follow. We use equation @ and convert the expression to Fourier space. Consider first only the h zz 
term, whose contribution to a\ we denote o\ a . The polarization gives (ef z ) 2 + (e zz ) 2 = sin 4 #fc, where 

k = (k,9k,(j>k) m spherical coordinates. Performing the angular k integrations then yields 
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a P,a 



432vr 
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TO 



TS 



TO 



dk k(n - t s )(t 2 - t s ) 



ji(fcn) h(kr 2 ) h[k(n - r 2 )] 
kr\ kr 2 [k(ri - r 2 )] 3 



The jz\k(r\ — T2)]/[^( r i — T 2)] 3 term represents a further suppression of short wavelength modes due 
to phase cancellations among different waves in the assumed isotropic stochastic background. Letting 
s = kr\ and q = T2/V1, we can simplify this expression to a double integral, 



a 



P,a 



864vr . 

^2-^T 



'/Ml [ ^TQ 



rs ~ rsrs + r s 



— + ln(q/q s ) 
2q 



W(q)dq 



(13) 



where qs = ts/tq and with u = (1 



W(q) 



s we define 

h(s) ji(qs) j 3 (u) 



s ds 



qs 



(14) 



Similarly, the contribution of the h z i terms of Eq. (0) is 



2 



1447T A 
~~J2 T 



G(q)dq 



(15) 



where 



G(q) 







2 1 

+ 5 


~ji(qs)~ 
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Ji(s)ji(qs) 
s qs 



3i(u) J2{u) 



(16) 



Integrating over angles gives a zero cross term, and so a 2 g = a 2 3 a + <x| 6 . Numerically, we find that 
ag = 5 x (zs = 1), 9 x 10 -6 (zs = 3). This is much smaller than the estimates for LSS |1C], 



6 x 10" 4 (z s = 1), 7 x 10" 4 (zs = 3). 

To estimate the relative deflection between rays at 9 = 0, we choose two directions (labeled A and B) 
separated at the observer by an infinitesimal angle 7, and find the rms difference between the deflections 

2\ V2 

We cannot evaluate this with 



due to GW's in these two directions, = y Pa(ts) — Pb(ts) 

the method used for LSS, which assumes that horizon size modes are negligible H. Instead we must 
calculate a\ Q explicitly and keep all the terms to lowest order in 7, i.e., quadratic order. These include 
terms which come from multiplying polarization components for the different directions A and B. The 
final result is 

= 4(tL + o\a + a\a b + a\a , (17) 



where 
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sds , h{q) = 

2 h(s)ji(qs) 
s qs 



h(q)dq , 



(18) 



ji(s)ji(qs) j4(n) s 3 ^ 




Then aA/3/7 = 7 x 10 -6 (z s = 1), 1.3 x 1(T 5 (z s = 3). By contrast, LSS gives a a^ 3 /l = 0-07 (z s = 
1), 0.14 (zs = 3). For LSS, the relative deflection is greatly increased by coherent deflections for short 
wavelength modes, but for GW's the effect of these modes is cut off by the redshifting as well as the 
temporal oscillations. We also used the relation krs^f -C 1 in the calculation of oa/3- The reason we 
find a oa/3 of order 7073 is that long wavelength modes overlap over the two light rays, and the relative 
deflection is small compared to the total deflection. Indeed, a Taylor expansion suggests that in general 
^a^/t ~ krsC/3, and kr$ ~ 1 is dominant for this GW spectrum. As shown above, the shear tensor 
(which is also used in weak lensing) is closely related to a A/3/ '7, and so the mean square ellipticity at a 
point induced by GW's is of order 10 -5 , again negligible compared to the few per cent expected from 
LSS (e.g., §). ' 

We can also try to derive general limits on GW's at astrophysical wavelengths from the induced shear. 
To obtain a limit on we compute the erA/3/7 produced by an isotropic background of GW's at a 
single wavenumber k. Note that for modes at a given k, we can use equations @ and (|j) even for 
short wavelengths (with a normalization factor, separate for each k), for times r long after matter- 
radiation equality. Since GW's at Horizon wavelengths are already strongly constrained by the CMB as 
noted above, we restrict our calculation to the case kr$ 3> 1, in which case the h z i terms in Eq. (pi) can 
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be neglected. We can estimate from Eq. @ that in order of magnitude <r| should equal At/(A;to) 4 , and 
thus that ~ At/(A;to) 2 . However, we find from the exact calculation that there is no term this 

large, only higher order terms in 1/(A;tq). We outline in the Appendix a mathematical argument showing 
this cancellation at small wavelengths. This result requires both the phase cancellations that come in 
averaging over an isotropic background, and also the oscillation with time of the GW's. With different 
assumptions, e.g., if we analyzed GW's from a particular source, which are then not isotropic, stronger 
limits may be possible. 

IV Proper motions induced by LSS and GW's 

We now consider the fluctuation of the angular deflection of image rays with time, and the resulting 
proper motion. If the deflection of image rays induced by LSS or GW's changes significantly during an 
observation of a lens system, then the slow shift in alignment between the lens and the source will change 
the impact parameter at the lens of a given ray from the source. The images will therefore move, and 
even tiny motions may be detected since the source motion is magnified if it is lensed by a primary lens. 
We first show that this effect is still expected to be too small to measure for LSS and for the GW power 
spectrum that we have considered above. However, given the weakness of existing limits on GW's at 
astrophysical wavelengths (Sec. I), we also consider possible limits on a general GW spectrum. 

Again we consider a single light ray from the observer out to some distance r$, in the absence of a 
primary lens (we consider the effect of a lens below). Given a ray with a fixed direction at the observer, 
its position x l (rs) at r$ moves with time, and it is this motion which we evaluate. In practice, we are 
interested in a fixed source at rs, in which case its apparent position will drift with the same speed but 
in the opposite direction. For LSS we have (e.g., fTo[) 

rrs 

x\r s ) = -2 / (r s - r)V^(r = r - r)dr, (19) 

Jr=0 

in terms of the Newtonian potential (or scalar metric perturbation) (j). We are now using the parameter r 
rather than r, since as time changes all comoving distances remain fixed. The only change is the time of 
evaluation of <j), and so to find dx 1 (r s) / dTQ from x l (rs) we simply replace 4>(t = tq — r) by <P(t = tq — r), 
with the partial time derivative in (j) taken at a fixed position. The rms value of dx l (rs)/dTQ depends on 
the power spectrum of <j), a quantity which has been estimated by various authors in connection with the 
Rees-Sciama effect on the CMB (e.g., [|^]). While the integrated deflection is dominated by short (~ 1 
Mpc) wavelengths, the LSS potential only evolves on a cosmological time scale. In an Einstein-deSitter 
universe, (j) is time independent in the linear regime of small density perturbations, but in this case too 
(j) becomes nonzero when nonlinear structure forms. In general, therefore, the proper motion induced by 
LSS is of order ct^/tq ~ 10~ 8 arcsec per year. For the gravity wave spectrum considered above, horizon 
size modes are dominant, and so here too the induced proper motion is of order ap/rQ, with a erg smaller 
by ~ 10 2 than for LSS. Any observed proper motion will thus be dominated by peculiar velocities of 
hundreds of km/s generated, e.g., by the velocity dispersion of stars in a galaxy or galaxies in a galaxy 
group or cluster. 

We now estimate the lensing limit on stochastic GW's in general, at any wave number k. VLBI 
observations can directly measure or limit proper motions, and this then implies a limit on GW's. Again 
we restrict ourselves to wavelengths with krs 3> 1, and consider first the apparent motion of a source 
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that is not lensed by a primary lens. The apparent motion due to GW's of a fixed object at distance rs 
is — d(5{rs) / d,TQ. Up to corrections of order 1/krs, the mean square of this motion is 



18vrA T 
5Pr n 4 



1 + - cos(2A;ro; 



(20) 



However, when there are both a lens and a source, a GW background will produce correlated proper 
motions in both, and the relative motion may be small. Limits from VLBI on proper motions in gravi- 
tational lenses were recently considered in Ref. p3] 1, and we proceed similarly. We may hope for strong 
limits because, in the presence of lensing, a proper motion of the source relative to the lens is magnified 
into a larger proper motion of the images. Furthermore, only a relative motion between images needs to 
be detected, as opposed to a more difficult measurement of motion with respect to an external reference 
frame, since if the source moves (relative to the lens), the different images do not all move together. In 
general, different values of the magnification matrix at the different image positions will produce relative 
motions between images of the same order of magnitude as the absolute motions. Moreover, pairs of 
highly magnified images generally have antiparallel motions [23]. 

To analyze how proper motion due to GW's may be magnified, we start from Eq. (|l2"|), and consider 



the same equation a time At later, when the deflections from GW's have changed. E.g., Oq L has changed 
to CiQ L + A^ 1 ), and a total change A in the observed 6' has been induced. Expanding the lens equation 
to first order in the small changes and solving for A, we obtain 



A,; = A, 



(2) 



A (2,3) 



A (2) 



(21) 



where the magnification matrix m( equals the inverse of b\ — dia{ ens (and is evaluated at 6' + otQL ~ ®ol)- 
Consider first the magnified term, A^ 2,3 ^ — A^ 2 \ Averaging over directions of A^ 2,3 ) — A^ 2 ) we obtain a 
result analogous to Eq. (^) for the mean square. Since M tJ is symmetric for a thin lens [O], it has two 
real eigenvalues m a and mf, (where the magnification M = \m a mb\). Letting 



M 



(to 2 + mf) 



1/2 



(22) 



we find that 



rms 



Ml 



A (2,3) 



A (2) 



M x rms 



£(2,3) _ £(2) 



(23) 



In Eq. 



3 3 

we may evaluate the rms on the right-hand side using a straight path (as in Sec. III). Letting 



CLTq CLTq dTQ 



(24) 



we find that in ^/? 2 s^> there is no term of order / (Iz^Tq) [as in Eq. (20)], but only higher order terms 
in 1/(/cto). Once again this small wavelength cutoff results from combining the time oscillation of GW's 
and the phase cancellations in averaging over an isotropic background (see the Appendix) , and as a result 
there is only a very weak limit on 
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V Conclusions 



Gravitational lensing is affected by perturbations to the homogeneous and isotropic background metric. 
Such perturbations, whether they are caused by LSS or GW's, may produce a number of effects on light 
propagation. One such effect is an overall shift in the angular positions of nearby objects, which is not 
observable. Another is a relative difference between the induced shifts in nearby light rays. This relative 
deflection manifests itself as a shear which may cause weak lensing and also affect strong lensing. A third 
effect is a fluctuation of the angular position of distant objects with time, leading to a directly observable 
proper motion. 

The actual amplitude of long wavelength modes of LSS and GW's is limited by the quadrupole 
anisotropy of the CMB. Even if both make comparable contributions to the anisotropy, LSS is dominant 
in its effects on lensing. This results from cancellations due to the time oscillation of short wavelength 
gravity waves, as well as the redshifting of their amplitude. For LSS, on the other hand, the effect of 
small coherent structures is amplified as the deflection executes a random walk. We find that the rela- 
tive deflection due to GW's is four orders of magnitude smaller than that of LSS, and is therefore not 
observable. 

The induced proper motions expected for LSS or for GW's generated in inflation are small compared to 
typical peculiar velocities, and thus are not observable. The motions are also too small to yield interesting 
limits on the energy density of GW's at shorter wavelengths. 

After this paper was submitted for publication, the bending of light by gravity waves was analyzed 
differently in Ref. [p4| ], for the case of short (subhorizon) wavelengths, in a nonexpanding flat space (i.e., 
neglecting the redshifting of the amplitude of GW's). That simplified analysis shows that the relative 
proper motion between two sources is small not only if they are at different redshifts along the same line 
of sight (in agreement with our calculation of {pL,s) in ^ ec - -^0' but also if they are separated on the sky 
by a small angle. The treatment presented in Ref. p4| changes quantitatively if expansion is included, 
but not qualitatively for GW's with a period short compared to the redshifting time scale (i.e., a Hubble 
time). 
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A Appendix 

In trying to set limits on GW's at short wavelengths kr$ S> 1, we twice encountered a weaker limit 
than simple dimensional analysis would suggest: Once in calculating ctap/j or shear in Sec. Ill, and 
then in estimating the magnified proper motion in Sec. IV. In this Appendix we outline the first of these 
calculations and show how this result emerges. The second calculation can be done similarly. 
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From Eq. (^) in the limit of short wavelengths (compared with the present Horizon), we derive 



2 /J_ 172 i ff ^_,.4 r o ^_/, „.\t^ ^cos(frri) cos(fcr 2 ) j 4 [A:(ti - t 2 )] 

s 



=^T ATk L s dTl j TS d ^-s)(r2-r s ){r,-n){r,-r 2 ) _ 



The ji[k{ji — T2)]/[k(ri — T2)] term comes from the angular k integrations, including the angular depen- 
dence of the polarizations and assuming an isotropic background. Letting x = kr\ and u = k(ri — t 2 ) 
(also xq = krQ, etc.) leads to 

864vr [ x ° (x-x s )(x -x) f x o-x ( x + M _ Xs )( 2 ; _ X _ U ) J4 ( M ) _ 

9 At / ax 5 / cm ; ^ -j— cos u(l + cos 2x) — sin «sm 2x\ 

k 2 r z s Jx s x 2 Jxs-x {x + u) z u 4 

We evaluate only the first cos u term here, since the other terms in the square brackets can be evaluated 
similarly. Note that dimensional analysis suggests that the x and u integrals should give a term of order 
1 (not larger, because of the oscillating integrand). 

To do the u integral we separate the smooth and oscillating parts and then repeatedly integrate by 
parts: We let w^ n \u) be the n-th indefinite integral of [cos u] j'4 (it)/u 4 with respect to u, and v^ n \u) the 
n-th derivative of (x + u — xs){xq — x — u)/(x + u) 2 with respect to u. For each n such that w^ n \u) 
converges for u — > ±00, we fix the arbitrary constant by w^{oo) + w^(— 00) = (Any constant will do 
for other n). Then the u integral equals a series of terms evaluated at the two limits of integration, 



00 

EM 

n=0 



Since the two series of terms can be handled similarly, we evaluate here only the u = xq — x terms. We 
do the x integration in the same way as the u integration. Thus we continue to integrate w^ n+l \x) with 
respect to x, and let v (n,m \x) be the m-th derivative of (x — xs)(xq — x)v^(xq — x)/x 2 with respect to 
x. The contribution to from the terms we have kept is 



864vr 



At (-l) n+1 {^ [n+m+2 ko)w [n ' m k^o)-w [n+m+2] (^o-x 5 )w [n ' m] (^)} • 



xl 

S n,m=0 

'(x) at x > 1 is of order x~^ n+m \ u>^(0) is or a constant, and we find that w^(x) at 
\x\ S> 1 is of order |x| n ~ 5 . This last fact, that w^(±oo) converges for the first few n, depends on the 
specific function w^(x) which in turn is determined by the two physical assumptions of time oscillation 
and angular averaging. The only term from the final sum that could give a contribution of order At jx 2 s 
is the n = m = term. We find that w^(0) is a nonzero constant, but since v^°'°\x) = identically, 
there is no term of this lowest order. 
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